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Abstract 

We discuss a problem concerning the ultrahigh energy neutrino propa- 
gation through the Earth. We present calculation of the neutrino-nucleon 
cross section at high energies, based on the unifed evolution equation at 
small X. We also show the solution of the transport equation for different 
neutrino fluxes originating from active galactic nuclei, gamma ray bursts 
and top-down model. 

1 Introduction 

The ultrahigh energy neutrino physics is an interdisciplinary subject which can 
address many vital problems in astrophysics, particle physics and geophysics, 
for a review see Neutrinos can traverse large distances without being 

disturbed and thus give us important information about distant astronomical 
objects like active galactic nuclei or gamma ray bursts. The ultrahigh energy 
neutrinos can have energies up to 10^^ GeV, much larger than currently accessi- 
ble at present colliders, like HERA. At these ultrahigh energies the structure of 
nucleon is probed at very small values of Bjorken x. It means that we are pos- 
sibly entering a region in which partons have large density. Therefore detailed 
knowledge of parton distribution function at very small x is vital to estimate 
the neutrino-nucleon cross section. This cross section rises strongly with energy 
and therefore neutrino propagation through the Earth can be affected by the in- 
creased interaction with matter at these ultrahigh energies. Thus it is necessary 
to consider the effect of attenuation of high energy neutrinos while traversing 
the Earth. This phenomenon could be used as a possible method of the "Earth 
tomography" by neutrinos. All these topics will be studied at large neutrino 
telescopes like AMANDA, NESTOR and ANTARES (for review see [§). 

In this paper we present a method of calculation the neutrino-nucleon cross 
section using gluon and quark distributions obtained from a unifed evolution 
equation for small x This equation embodies both DGLAP and BFKL evo- 
lutions on equal footing as well as important subleading effects via consistency 
constraint This constraint limits the virtualities of the exchanged gluon 
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momenta in the ladder to their transverse components. We then extrapolate 
the results to ultrahigh energies and compare the predictions with the other, 
based on standard global parton analysis 0,1^. 

We use the resulting cross section as an input to the transport equation [Q 
for the neutrino flux penetrating the Earth. We solve this equation for different 
incident angles as well as different input neutrino fluxes originating from active 
galactic nuclei, gamma ray bursts and a sample top-down model. 

The details of the calculation has been presented elsewhere M . 



2 Sources of ultrahigh energy neutrinos 

Neutrinos have the advantage that they are weakly interacting with matter and 
they are hardly absorbed when travelling along large distances. On the other 
hand cosmic rays are being absorbed through the following interactions: 

• pair production: jcR. + Ibr e+e~ 

• inverse Compton scattering: ccb. + Ibb. e + 7 

• photoproduction of pion: pcR + "fBR — > p + Nt: 

• nuclei fragmentation by photo-pion interactions 

These reactions cause that the spectrum of the cosmic rays has the GZK cutoff 
around lO^^eV. Therefore neutrinos are the best candidates for supplying 
information about distant objects. The possible sources of highly energetic 
neutrinos arc: 

• active galactic nuclei 

• gamma ray bursts 

• top-down models 

Active galactic nuclei are the most powerful sources of radiation in the whole 
Universe. The engine of the AGN is the supermassive black hole with mass 
bilion times larger than the mass of the sun. Surrounding the black hole is the 
accretion disc usually accompanied by two jets. The spectrum of the emitted 
photons spreads from radio waves to TeV energies. The jets are the sources 
of the most energetic gamma rays, since the particles (electrons and possibly 
protons) are accelerated in blobs along the jets with Lorenz factor 7 10. 
Electrons loose energy via synchrotron radiation thus producing very energetic 
photons. If protons are also accelerated then they can interact with the am- 
bient photons producing pions. Consequently strong flux of neutrinos will be 
produced 0, |ll|, 
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Gamma ray bursts are also very probable sources of neutrinos. Although 
the underlying event of GRB's is not entirely known the present knowledge is 
consistent that the bursts are produced as a relativistically expanding fireball 
which intitial radius was around 100km. The original state which is opaque to 
light expands in a relativistic shock with 7 ~ 300 to the point where it becomes 
optically thin and produces intensive gamma ray spectrum. As in the case of 
AGN, the acceleration of protons can result in the production of neutrinos be- 
cause they will photoproduce pions and in the end neutrinos will emerge [p^ . 

Top-down models are the most speculative scenarios of producing neutrino 
fluxes at ultrahigh energies. In these models one assumes that the particles are 
not accelerated but are rather produced as a result of the decay of the supermas- 
sive particles ^ 10^^ — lO^^GeV. These particles could emerge as a decay of 
some topological defects: superconducting strings or magnetic monopoles. The 
top-down models produce quite hard spectrum of neutrinos extending beyond 
lO^^GeV and they were proposed as a possible solution to the GZK cutoff ]T^ . 



3 Neutrino - nucleon cross section 

The dominant interaction for neutrinos is the neutrino-nucleon interaction. It 
has the largest value of the cross section, and is dominating over the interaction 
with electrons. There is however one exception: resonant W production in i7e e~ 
interaction at E^, = 6.4 x 10^ GeV. At this energy this process domintates by 2 
orders of magnitude over other contributions. The neutrino - nucleon interaction 
can be visualised in Fig. |l|. 
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Figure 1: Deep inelastic scattering. 
Here / is the incoming neutrino of four momentum k and I' is the outgoing 
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neutrino or charged lepton with four momentum k' . N is the target nucleon with 
four momentum p and X is the arbitrary hadronic final state. This deep inelastic 
scattering process can be described using two standard variables: q"^ = —Q"^ < 
is the four momentum squared of the exchanged vector boson, and x = is 
the standard Bjorken scaling variable. At high energies (up to lO^^GeV) , very 
small values of x can be probed because, 

x ^ .^hL_ ^ 10-8 (1) 

where My is the mass of the heavy vector boson, is the nucleon mass 
and is the energy of the exchanged vector boson. These are values of x which 
are not accessible at present electron-proton colliders (for example HERA goes 
down to 10^^ at fairly low Q^). Therefore a detailed knowledge of parton 
distributions at these small values of x are required. We propose to use the 
unifed BFKL/DGLAP evolution equation Q]. We shall start at first with the 
pure leading order BFKL equation for the unintegrated gluon distribution 
function f{x,k'^) in the foUowng form: 

fix e) f^^ix e)+a,k^ rdfc-^ f fix/z,k'^)~fix/z,e) f{x/z,k^) 

(2) 

where as = NcOis /t: and k — kT,k' — kip denote the transverse momenta of the 
gluons, see Fig. 2. The term in the integrand containing f{x/z, k'^) corresponds 
to real gluon emission, whereas the terms involving f(x/z,k'^) represent the 
virtual contributions and lead to the Reggeization of the i-channel exchanged 
gluons. The inhomogcncous driving term Z*^*^^ is the input function will be 
specified later. This is the leading order in ln(l/a;) equation which gives the very 
well known intercept for the gluon A — In 2. The next to leading contribution 
to the t-channel exchange at high energies has been already calculated 
It yields however very large correction to the intercept making it physically 
unreliable. It occurred, that the resummation of the subleading effects should 
be performed in order to get physically reliable results |l^ . Different forms 
of resummation has been already proposed in the literature. It appears that 
the imposition of the consistency constraint |^ which limits the phase-space 
available for the real emission term provides with the partial resummation of 
the subleading effects. This constraint arises from the fact that in the high 
energy limit the virtualities of the exchanged momenta are dominated by their 
transverse parts. By imposition of this constraint one obtains nice physical 
picture in which the subleading effects are resummed by the limitation of the 
phase space. The result for the gluon intercept yields reasonable value which 
is stable, i.e. does not become negative. Second improvement to the equation 
(H) is the inclusion of the DGLAP terms which are important for large values 
of X and the overall normalisation of the resulting gluon distribution function. 
We do also include the quark driving term in eq. 0. The resulting evolution 
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equation for the gluon has the following form: 




We specify the driving term in the following form 

Let us note that the inhomogenious term has been entirely specifed in terms 
of the standard integrated gluon distribution function at the infrared cutoff /cq 
and that the BFKL/DGLAP evolution only takes place above this cutoff. The 
last term is the contribution of the quark distribution to the gluon evolution, 

E = ^ x{q + q) = Y,{S, + V,) (5) 
9 9 

where S and V denote the sea and valence quark momentum distributions. The 
gluon, in turn, helps to drive the sea quark distribution through the g qq 
transition. Thus equation has to be solved simultaneously with an equivalent 
equation for T,{x, k^). We use the kr factorisation formula |T^, see Fig. ^ as a 
basis for our evolution equation for the quark distribution. 

S,{x,Q') - j ^Sl^{z,k^Q')f[^^,k') (6) 

where S^°^ describes the quark box (and crossed-box) contribution shown in 
Fig. 2 and can be interpreted as a partonic structure function. It has been 
shown (see for example that the ln(l/a:) effects are also resummed and 

play important role in the kx factorisation prescription. In order to get the 
complete set of evolution equations we also have to add quark-quark splittings 
(which are small numerically anyway) and the valence quarks, which we take 
from the set of parametrisations. Thus our complete equation for the singlet 
quark distribution reads as follows, 

E(X, k') = 5„on-p(x) + E /" T ^9™(^' = 0' ^')? -9 ^o) 

q ^ 
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Figure 2: Diagrammatic representation of the /cT-factorization formula. At 
lowest order in as the gauge boson-gluon fusion processes, Vg —>■ qq, are given 
by the quark box shown (together with the crossed box). The variables k, k 
and k' denote the transverse momenta of the indicated virtual particles. 

(7) 

where a — {1 + Am^/Q'^)^^ and V — x{uy + d„). Equations (^) and (0) form a 
set of coupled integral equations for the unkown functions f{x, k"^) and S(a;, k"^). 
We solve them assuming simple parametric form of the inputs: 

xg(x,kl) = N{l-xf 
^non-p(x) = C7p(l - (8) 
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We have used this set of coupled equations to calculate F2 structure function 
at HERA, and we have therefore fixed the values of the free parameters. We 
then used the resulting parton distribution functions in order to calculate the 
neutrino - nucleon cross section at high energies. We calculate the cross sections 
from the usual formula: 

d^a--^ GfME ( Mf V + 

2 ^2 [x, Q ) 

(9) 



dxdy TT V + M 



^i^i;(x,Q2) + y(l-|) xF,^{x,Q') 



where Gp is the Fermi coupling constant, M is the proton mass, E is the 
laboratory energy of the neutrino and y — Q'^ /xs. The mass Mi is either Mw 
or Mz according to whether we are calculating charged current (CC) or neutral 
current (NC) neutrino interactions. Functions F^, F£ and xF^ are of course 
usual structure functions which can be calculated from the parton distribution 
functions. 

In Fig. 3 we show the plot of a{vN) as a function of energy of the neutrino E^. 
One observes strong rise, nearly 8 decades with increasing energy from 10 to 
lO^^GeV. The shape of the curves up to lO^GeV is determined by the valence 
quarks whereas beyond lO^GeV everything is driven by the sea quarks. One 
also notes that the charged current contribution is dominating over the neutral 
current by factor 3. 



4 Transport equation 

Neutrinos at ultrahigh energies can be quite strongly attenuated when traversing 
the Earth. Apart from standard absorption neutrinos can undergo regeneration 
due to neutral current interactions. Charged current interactions remove neu- 
trinos from the flux, but neutral current interaction cause neutrinos to reappear 
at lower energies. These both effects can be calculated using the transport 
equation proposed by : 

dI{E,T) ... . , f dy danc[E',y) 



dT 



J ^-y dy 



where cttot = ""cc + "'nc ^-nd where y is, as usual, the fractional energy loss 
such that 

E' = j^. (11) 

1 - y 

The variable t is the number density of nucleons n integrated along a path of 
length z through the Earth 

t(z) = / dz'n{z'). (12) 
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Figure 3: The total vN cross section and its decomposition into contributions 
from charged current and neutral current plotted as a function of the laboratory 
neutrino energy. 

The number density n{z) is defined as n{z) — Na p{z) where p{z) is the density 
of Earth along the neutrino path length z and Na is the Avogadro number. 
The number of nucleons r encountered along the path z depends upon the 
nadir angle 9 between the normal to the Earth's surface (passing through the 
detector) and the direction of the neutrino beam incident on the detector. 

In order to calculate the change of the intensity of the flux with the incident 
angle one needs to know the density profile of the Earth. We have used the 
model by A. Dziewohski Using this parametrisation of the Earth density 
and the cross sections calculated from the unifed BFKL/DGLAP equation. In 



8 



Regenemtion+absorption - 

Absorption 

Atmospheric 

9=0° \ - 


\ AGN (Protheroe) 

6=0° '\\ 


— 1 1 1 Mil 1 1 1 Mil 1 iTmi^ 1 1 Mill 1 1 Mil 

6=40° 


— 1 1 1 Mil — 1 1 1 Mil — rr?ftm] — 1 mill — 1 1 Mill 
6=40° 


— 1 1 1 Mil 1 1 1 Mil 1 1 Mill 1 1 Mil 


— 1 1 1 Mil — 1 1 1 Mil — 1 1 1 Mil rrTrrml — i i iiiii 


6=80° - 


6=80° ~^^^^^^^^^^^^~^^N^ 



34567345678 

10 10 10 10 10 10 10 10 10 10 10 

Energy [GeV] 



Figure 4: The shadowing factor S of (|T^) for two different initial neutrino 
fluxes incident at three different nadir angles on a detector. The angle 9 — Q° 
corresponds to penetration right through the Earth's diameter. The two curves 
on each plot show the shadowing factor with and without NC regeneration 
included. 



Fig. 4 we show the shadowing factor, 

S{E.r) = ^ (13) 

where I^{E) = I{E, r) is the initial flux at the surface of the Earth. We present 
iS* as a function of the energy and for different incident angles. The curves 
exhibit strong suppression for large paths in matter and the energies above 
lO^GeV. Also, the curves corresponding to the flux from AGN jl^, show that 
the regeneration is important for flat fluxes and large paths. The same effect is 



nearly negligible in the case of steeply falling atmospheric spectrum 22 
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Figure 5: The initial flux Iq{E) and the flux at the detector I{E) for three 
different nadir angles corresponding to three models for AGN neutrinos |jl^, 
|l2|. The background atmospheric neutrino flux is also shown. All the fluxes 
are given for muon neutrinos. The corresponding fluxes from | pO| , p^ , 12 1 were 
given originally for muon neutrinos and anti- neutrinos, and their value has been 
divided by factor 2. 



In Figs. 5 and 6 we show complete simulation for different incident fluxes 
corresponding to the active galactic nuclei jl^, |ll], |l^, gamma ray bursts |Q 
and a sample top-down model ||lj]. Similar effects of strong attenuation are 
observed for large energies ~ lO^GeV^ and large paths. This suggests that one 
will have to choose suitable angle of the observation in order to avoid large muon 
background from the atmosphere and yet be able to detect the highly energetic 
neutrinos which arc likely to be absorbed by the matter in Earth. 
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Figure 6: The initial flux Iq{E) and the flux at the detector I{E) for three 
different nadir angles corresponding to the model of gamma ray burst and 
the top-down model p^ . The background atmospheric neutrino flux is also 
shown. All the fluxes are given for muon neutrinos. 

Summary 

In this paper we have examined the interactions of ultrahigh energy neutri- 
nos with matter. As the cross section rises with increasing energy the absorption 
by matter becomes interestingly large. We have seen that the Earth becomes 
essentially opaque to ultrahigh energy neutrinos. We have calculated neutrino- 
nucleon cross section using the unifed DGLAP /BFKL evolution equations which 
treat leading \n[Q'^) and ln(l/a;) on equal footing. These equations also resum 
important subleading effects in ln(l/a;) via imposition of consistency constraint. 
We believe that this form of the evolution is most apropriate in the regime 
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where the parton density is large. The results for the cross section has been 
compared with the other based on the standard global parton analysis. In that 
way the uncertainty due to parton density extrapolation has been diminished to 
40%. We have then used the resulting cross sections and solved the transport 
equation for the neutrinos travelling through the Earth. We have found that 
the attenuation is large for high energies , above lO^GeV and large paths in 
matter. We have also found that the regeneration due to neutral current inter- 
actions becomes important for flat spectra (like these originating from active 
galactic nuclei) and large paths. We have simulated the penetration through 
the Earth for different neutrino fluxes: active galactic nuclei, gamma ray bursts 
and top-down models. The large attenuation effect reduces substantially the 
flux for small nadir angles. There is however a window for observation of AGN 
fluxes by the km"^ detectors. 
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